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Abstract 

Universality in isotropic, abelian and non-abelian, sandpile models is ex- 
amined using extensive numerical simulations. To characterize the critical 
behavior we employ an extended set of critical exponents, geometric features 
of the avalanches, as well as scaling functions describing the time evolution of 
average quantities such as the area and size during the avalanche. Comparing 
between the abelian Bak-Tang-Wiesenfeld model [P. Bak, C. Tang and K. 
Wiensenfeld, Phys. Rev. Lett. 59, 381 (1987)], and the non-abelian models 
introduced by Manna [S. S. Manna, J. Phys. A. 24, L363 (1991)] and Zhang 
[Y. C. Zhang, Phys. Rev. Lett. 63, 470 (1989)] we find strong indications 
that each one of these models belongs to a distinct universality class. 
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I. INTRODUCTION 



Sandpile models have been studied extensively during the last decade as a paradigm 
of self-organized criticality (SOC). This concept, introduced by Bak, Tang and Wiesenfeld 
(BTW) []I],0]> stimulated numerous theoretical [[3-14], numerical fl5|-22| and experimental 
studies P3|-f28|[. In sandpile models, which are defined on a lattice, grains are deposited ran- 
domly until the height at some site exceeds a threshold, and becomes unstable. Grains are 
then distributed to the nearest neighbors. As a result of this relaxation process neighboring 
sites may become unstable, resulting in a cascade of relaxations called an avalanche. It was 
observed that these models are self-driven into a critical state which is characterized by a 
set of exponents |l]||. This set includes distribution exponents that describe the distribu- 
tion of quantities such as avalanche size, area, and lifetime, and the geometric exponents 
which relate various properties of the dynamics ||21|| . Additional sandpile models which dif- 
fer from the BTW model in the dynamic rules were introduced and studied. These include 



the Manna model [18, 20 1 , in which the dynamics during an avalanche is stochastic and the 



Zhang model fll5| , in which the dynamic variable is continuous. In order to understand the 
basic mechanism of SOC, numerous attempts were made to assign the various models to uni- 
versality classes. To this end, Vespignani, Zapperi and Pietronero introduced the fixed scale 
transformation which is a real space renormalization group (RG) approach P0| , |3~T| . They 
applied this approach to the BTW and Manna models and concluded that these two models 
belong to the same universality class. Later, a comparison between the critical behavior of 
BTW and the Manna models was performed, using extensive numerical simulations and an 
extended set of exponents |32fl . The results provide strong evidence that the two models 
belong to different universality classes. 

Recently, Diaz-Guilera and Corral applied the dynamic RG technique to study critical 
behavior in sandpile models. Their conclusion was that the BTW and the Zhang models be- 



long to the same universality class p3|-|35|] . A similar conclusion was reached by Liibeck |36 
based on numerical simulations. 
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In this paper we apply an extended set of tools for the characterization of critical behavior 
in sandpile models and for their classification into universality classes. These tools include 
measures that characterize avalanches as a whole such as the distribution exponents, the 
geometric exponents and the avalanche structure. We also introduce a new measure, based 
on scaling functions, describing the time evolution of the area, size and energy during an 
avalanche. Combining all these tools, we find strong indications that the BTW, Manna and 
Zhang models belong to three different universality classes. 

The paper is organized as follows. The models are introduced in Sec. II and the measures 
for their classification are presented in Sec. III. The simulations and results are given in Sec. 
IV, followed by a summary and conclusions in Sec. V. 

II. THE MODELS 

Sandpile models are defined on a d dimensional lattice of linear size L. Each site i is 
assigned a dynamic variable E(i) which represents some physical quantity such as energy, 
stress etc. In a critical height model a configuration {E(i)} is called stable if for all sites 
E(i) < E c , where E c is a threshold value. The evolution between stable configurations is by 
the following rules: 

(i) Adding energy. Given a stable configuration {-E(j)} we select a site i at random and 
increase E(i) by some amount 5E. When an unstable configuration is reached rule (ii) is 
applied. 

(ii) Relaxation rule. If E(i) > E C) relaxation takes place and energy is distributed in the 
following way: 

£(i)-£(i)-£A£(e) 

e 

(1) 

E(i + e) ^£(i + e) + A£(e), 
where e are a set of (unit) vectors from the site i to some neighbors. As a result of the 
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relaxation the dynamic variable in one or more of the neighbors may exceed the threshold. 
The relaxation rule is then applied until a stable configuration is reached. The sequence of 
relaxations is an avalanche which propagates through the lattice. 

It was shown before, that the parameters 8E and E c are irrelevant to the scaling be- 
havior |^ J53|rj^ |. This indicates that the critical exponents depend only on the vector AE, 



to be termed relaxation vector. For a square lattice with relaxation to nearest neighbors 
(NN) it is of the form AE = (En, Ee, E$, Ew), where En, Ee, Es and Ew are the amounts 
transferred to the northern, eastern, southern and western NN's respectively. In the BTW 
model, E c = 4, 5E = 1 and AE = (1, 1, 1, 1). If an active site with E(i) > E c is toppled, it 
would not become empty after the topple had occurred. In the Zhang model |15[ , for which 



E c = 1 and < 5E < 1, the relaxation vector is given by (b, b, b, b), where b = _E(i)/4 and 
E(i) is the amount of energy in the active site before the topple had occurred . Obviously, 
the site i remains empty after toppling. 

In a random relaxation model J57J a set of neighbors is randomly chosen for relaxation. 
Such a model is specified by a set of relaxation vectors, each vector being assigned a prob- 
ability for it application. For example, a possible realization of a two-state model includes 
six relaxation vectors (1, 1, 0, 0), (1, 0, 1, 0), (1, 0, 0, 1), (0, 1, 1, 0), (0, 1, 0, 1) and 
(0, 0, 1, 1), each one applied with a probability of 1/6. A four state model would include 
relaxation vectors such as (4, 0, 0, 0), (3, 1, 0, 0), (2, 2, 0, 0), (2, 1, 1, 0), (1, 1, 1, 1), (0, 
2, 1, 1), etc., applied at different probabilities while maintaining the four fold symmetry of 
the relaxation rule. A time step (of unit time) is defined as the relaxation of all the sites 
having E(i) > E c , after the completion of the previous time step. A model is said to be 
abelian if the configuration after the avalanche, does not depend on the order in which the 
relaxation of the active sites was performed. The BTW model was shown to be abelian [17] . 



The Manna models |18|,23 are not abelian because they contain a random choice of the 
toppling direction. As a result, they develop different scenes of toppling that depend on 
the order of relaxation of the active sites in a single time step. The Zhang model is also 
non-abelian. This can be seen from the following example, when two active NN sites are 



toppled within the same time step. The site that was toppled last remains empty while the 
other one is non-empty. This shows that the final configuration depends on the order. 



III. MEASURES FOR CLASSIFICATION 

Avalanches have various properties which can be measured in a simulation: size, area, 
lifetime, linear size, and perimeter. The size (s) of an avalanche is the total number of 
relaxation events that occurred in the course of a single avalanche. The area (a) is the 
number of sites in the lattice where relaxation occurred. The lifetime (t) of an avalanche is 
the number of time steps (defined above) that took place during the avalanche. As for the 
linear size of an avalanche, one can use the radius of gyration (r) of the cluster of sites where 
relaxation occurred. A site inside the area a is defined as a perimeter site if it has a nearest 
neighbor where no relaxation took place. The perimeter (p) is the number of perimeter 
sites. Thus we have a set of variables {s,a, t,r, p} which characterize an avalanche. The 
avalanche variables have probability functions which are assumed to fall off with a power law 
defined by P(x) ~ x x ~ Tx , where x 6 {s, a, t, r,p} and the exponents r x are called distribution 
exponents. These variables also scale against each other in the form 

y~x lyx , (2) 

for x, y G {s,a,t,r, p} and the exponents 7 y:r are called geometric exponents. The exact 
definition of the geometric critical exponents ^ xy is in terms of conditional expectations 
values: ~ x lyx [21|. The exponents are not independent. It is shown in [^2[ that they 



satisfy the scaling relations 



lyx = 7™ (3) 



and 



(4) 
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Apart from looking at the critical exponents we have also examined the structural features 
of the avalanches. We define on the lattice a function that specifies the number of toppling 
events, /(i), for each site i during a single avalanche. In a two dimensional (2D) model this 
function takes the form of a hilly terrain, with discrete heights. The nth terrace is then 
defined as the set of sites for which /(i) > n. This height profile of an avalanche can be 
described by drawing the area of the nth terrace, A n vs. n. Furthermore, we can consider 
A n as an avalanche variable, just like the size or the area. This gives us new variables that 
can be measured and compared between the various models . 

The avalanche properties introduced so far, such as the area and size, characterize an 
avalanche as a whole, and are measured only after the avalanche is completed. Here we 
introduce measures to characterize the time evolution during an avalanche. Measured vs. 
time, for a single avalanche, these measures exhibit an irregular form. However, if we average 
them over a large number of avalanches, a typical shape emerges. We now introduce three 
such measures: 

1. The avalanche area a c (t), namely the number of sites where at least one relaxation 
occurred during the first t time steps of the avalanche. As the avalanche is completed 
a c (t) evolves to the area a of the avalanche. We also define a(t) as the time derivative 
of a c (t), according to a(t) = da c (t)/dt where da c (t)/dt = a c (t+ 1) — a c (t). The variable 
a(t) gives the number of sites that at time t became active for the first time (and are 
to be toppled in the next time step). As the avalanche evolves to an end we find that 
the avalanche area satisfies a = Y%=o x a (^)> where t max is the avalanche time, and the 
t — step consists of the deposition event that initiated the avalanche. 

2. The number of active sites s(t) (namely sites that are to topple in the next time step), 
as a function of time. As the avalanche evolves to an end we find that s = Y%=o x s (0 
is the avalanche size. 

3. Considering a sandpile in a gravitational field, we define the potential energy U(t) = 
X)iUi(i), where Ui(t) is the energy at site i given by : 
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fE(i) 

Ui(t) = / mgh(t)dh. (5) 



o 



The values of m and g are irrelevant and we take mg = 1. The potential energy defined 
here, turns out to be proportional to E(\) 2 , contrary to the ordinary definition of energy 
in sandpiles, which is linear in E(i). For the case of a discrete dynamic variable, a site 
i having energy E(i) will have a potential energy Ui(t) = E(i)(E(i) + l)/2. 

Averaging the functions a(t), s(t) and u{t) over a large number of avalanches we obtain 
the functions A(t), S(t) and U(t), respectively. According to the dynamic scaling assump- 
tion, each one of these functions can be written in the general scaling form : 

X(t) = K x (t) x a -f x (J^j (6) 

where 

(t) - ZttXW 

{t)x ~ E^m (7) 

and X E {U, 5, A}. The scaling function /x(/-0? where /x = tj {t)x satisfies the sum rules 

roc poo 

/ fx(v)dfJ. = / nfx{n)dn = 1. (8) 

JO JO 

The shape of the scaling function, and the values of the exponents ax can be used to distin- 
guish between universality classes. Moreover, the relation between these scaling functions 
can be used as a further tool. For example, if fx(n) and /y(/i) coincide for one model and 
are different in another model, it indicates that these two models do not belong to the same 
universality class. The dependence of (t) x on the system size L is given by 

(t)x-L^. (9) 



IV. SIMULATIONS AND RESULTS 



Having defined the three models and the measures used for their characterization, we 
now describe the computer simulations. We have used open boundary conditions and system 
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sizes up to 512 2 , with 10 6 to 10 8 grains dropped, in two dimensions (2D). For each run we 
ascertained, before collecting data, that the dynamics has reached the critical state by 
applying Dhar's burning algorithm |7]], or by starting with a configuration belonging to the 
critical state. For each model we have calculated all the measures of classification mentioned 
in Sec. III. 

The exponent t s that describes the avalanche size distribution was measured for the 
BTW, Manna four-state and Zhang models (Fig. [l]) finding good agreement with previous 



results |32] , |36[| . The exponent t s for the three models, as a function of the inverse system 
size, is given in Fig. 0. These results do not provide a reliable extrapolation of t s to the 
infinite system limit. However, they strongly indicate that the curves converge to different 



values of t s as L — > oo. As t s exhibits relatively strong dependence on the system size [p2[ , it 
cannot be used as the primary tool for classification of models, but only to provide additional 
evidence. 

The geometric exponents "f xy are only weakly dependent on the system size, and turn out 



to be very useful for classification of sandpile models P2] . In our simulations we examined 
the exponents ■j sa , , Jst and j at . This was done by drawing on a log-log scale quantities such 
as the average avalanche size £"[s|a] for a given area a, where 7 sa is given by the slope of the 
straight line section (Fig. §). The values of the geometric exponents for the BTW and the 



Manna four-state models, are in agreement with previous simulations |32j . Our simulations 
of the Zhang model, for system size up to 512 2 , showed that these variables, 7 sa ,7 s t and 
7a£, are not scale invariant, (Fig. ^j. Although for small avalanches the scaling behavior for 
the Zhang model resembles the distribution of the same variables in the BTW model, for 
large avalanches, we find a bend and a sudden change in slope. This bend is seen for all the 
system sizes that were checked. For large avalanches, the geometric exponents are clearly 
different from the values recorded for the BTW model and for the Manna four-state models. 
This puts in question the previous assignment of the Zhang model to the universality class 
of the BTW model CTTOim 



To further characterize the avalanche structure we examined the function /(i), that 



provides the number of toppling events at site i during the avalanche (Fig. |j). For the 
BTW model, we observe a shell structure in which all sites which relaxed at least n + 1 
times form a connected cluster with no holes which is contained in the cluster of sites which 
relaxed at least n times [[7|,p9|. The Manna four-state model exhibits a random avalanche 
structure with many peaks and holes [J32| . In between we find the Zhang model, which shows 
an avalanche structure which is mostly shelled, but is different from the BTW picture by 
having several peaks and holes, but not as many as in the Manna four-state model. 

To obtain a more quantitative characterization of the terrace structure we chose typical 
large avalanches for each of the three models and plotted the terrace number n as a function 
of its area (Fig. |5|). The results for the Manna model are well fitted by an exponential, while 
a nearly linear decay is observed for the BTW and Zhang models. We have also measured, 
for the three models, the distribution exponents r a (n) for terraces no. 2, 3, 4 and 5 (Table |). 
The results show quantitative differences in the avalanche structures, between the three 
models. In all cases r a (n) decreases as n is increased. The differences are significant with 
the lowest r a (n) for Zhang, intermediate for BTW and highest for Manna models. To obtain 
a more complete characterization of critical behavior in sandpile models we also examine the 
time evolution of the energy, avalanche size and area during the avalanche, averaged over 
a large number of avalanches. Combining results for system sizes L = 128, 256 and 512 we 
draw the scaling functions fu( T ), fs( T ) an d /a( t ) that describe the averaged time evolution 
of the energy, number of active sites, and area growth during the avalanche. The scaling 
functions for the BTW, Manna and Zhang models are shown in Figs. ||, ^ and || respectively. 
For the BTW and Manna models we find excellent data collapse indicating scaling behavior. 
No such scaling is found for the Zhang model, indicating that it lacks some of the features 
of a critical system, which are found in the BTW and Manna models. For the BTW model 
(Figs |]) we observe that all the three scaling functions fu( T ), fs( T ) an d Ja{t) are identical, 
so the system is basically described by a single scaling function. For the Manna model, 
we find that fu{ T ) and /s(t) coincide, while /a(t) is different. For the Zhang model we 
find that there is no data collapse and there are no scaling scaling function. Interestingly, 
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for each system size the functions fuij) and fsij) are identical, while is different. 

The average times are found to depend on the system size according to (t)x ~ L /3x where 
X E {U, S, A}. For the BTW model (3 V = 1.51,/3 S = 1.43 and (3 A = 1.31; for the Manna 
model j3jj = 1.52,/?s = 1.53 and (3a = 1.48; for the Zhang model one can approximate these 
exponents by values (3u = 1.5, (3$ — 1-46 and Pa = 1-36 but there is a significant deviation 
from a straight line. 



V. SUMMARY AND CONCLUSIONS 

We have studied universality in isotropic, abelian and non-abelian, sandpile models us- 
ing a combination of extensive numerical simulations and an extended set of measures to 
characterize these models. In particular, we focused on the BTW model (which is abelian, 
deterministic and isotropic), Manna model (non-abelian, stochastic and isotropic on aver- 
age) and the Zhang model (non-abelian, deterministic and isotropic). For each model we 
have calculated the critical exponents which characterize an avalanche as a whole. These 
include the distribution exponents r x , that characterize the distribution of quantities such 
as avalanche size, area and lifetime, and the geometric exponents 7^, that relate the scaling 
properties of different quantities. The geometric exponents 7^ are particularly useful for 
classification due to their weak dependence on the system size. Comparing these exponents 
we find clear indications that the BTW and Manna models belong to different universality 



classes, in agreement with previous simulations |32fl . As for the Zhang model, the geometric 
critical exponents are not well defined. For all system size examined the functions 
vs. x, where x G {s,a,t}, from which the exponents ^ xy are obtained exhibit domains with 
different slopes for small and large avalanches. The small avalanche behavior is similar to 
the BTW results, while the large avalanche behavior is different from both the BTW and 
Manna models. 

The avalanche structures of the three models are found to be significantly different. The 
BTW avalanche structure is the most regular, the Manna structure is the most irregular 
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and the Zhang avalanche structure is intermediate. 

We have also examined measures of the dynamics during the avalanche. We found scaling 
functions for the time evolutions of the energy fu{ T ), number of active sites /s(t), and the 
rate of area growth /a(t) in the BTW and Manna models. For the BTW model, all three 
scaling functions coincide, while for the Manna model only the first two coincide. This is a 
qualitative different that further strengthen our conclusion that the two models belong to 
different universality classes. For the Zhang model these functions do not exhibit scaling 
behavior. This is a further indication that the Zhang model lacks some essential features 
of critical behavior, which appear in the BTW and Manna models, and thus belong to a 
different universality class. In fact, the only unambiguous scaling features of the Zhang 
model are given by the distribution exponents t x . We thus conclude that the BTW, Manna 
and Zhang models belong to three universality classes. 

Our results disagree with the conclusions of a number of recent studies. Liibeck stud- 
ied the scaling behavior in the BTW and Zhang models using extensive numerical simula- 
tions [5^]. Relying only on the distribution exponents he concluded that the BTW and Zhang 



models belong to the same universality class. As we demonstrated above, the distribution 
exponents provide very limited characterization of the scaling behavior. Therefore, these 
exponents alone do not provide enough information to support a conclusion that two mod- 
els belong to the same universality class. Moreover, the distribution exponents are strongly 
dependent on the system size. One should be careful in interpreting the results of the finite 
size analysis done in |36[], based on an assumed size dependence which is not substantiated 



theoretically. Corral and Diaz-Guilera derived nonlinear partial differential equations based 



on the microscopic evolution rules of the BTW and Zhang models ||35|| . Using a dynamic RG 
approach they studied and concluded that the two models belong to the same universality 
class. Vespignani, Zapperi and Pietronero used a real space RG approach and concluded 
that the BTW and Manna models belong to the same universality class [[JIJ. The failure 
of these approaches to distinguish between the universality classes indicates that some key 
ingredient of the dynamics is not taken into account. We speculate the these ingredient may 
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be related to symmetries such as the abelian symmetry, as well as properties such as the 
deterministic vs. stochastic nature of the avalanche dynamics. 
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TABLES 

TABLE I. The distribution exponent r a (n) for the areas of the terraces no. n = 2,3,4,5, in 
a 2D sandpile of size 128 2 . For all three models, r a (n) tends to decrease as the terrace order n 
increases. The differences between the models are significant, with the lowest exponents for Zhang, 
intermediate for BTW and highest for the Manna model. 
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Exponent model 

BTW Manna four-state Zhang 

r (2) 2.05 ± 0.03 2.16 ± 0.03 1.98 ± 0.03 

r (3) 2.00 ± 0.04 2.11 ± 0.04 1.85 ± 0.04 

r (4) 1.95 ± 0.05 2.09 ± 0.05 1.77 ± 0.05 

r (5) 1.91 ± 0.05 2.00 ± 0.05 1.74 ± 0.05 
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FIGURES 

FIG. 1. Avalanche size distributions in the BTW, Manna four-state and Zhang models. The 
exponents are r s = 2.09 ± 0.005 for BTW, t s = 2.23 ± 0.01 for Manna and r s = 2.25 ± 0.01 for the 
Zhang model. System size is 128 2 . 

FIG. 2. The exponent t s of the avalanche size distribution as a function of the inverse system 
size (1/L) in the BTW, Manna four-state and Zhang models. The lines are guides to the eye. 
Although it is hard to reliably extrapolate from these results to 1/L — > 0, this graph strongly 
indicates that the curves converge to three different values. 

FIG. 3. Geometric critical exponents for the BTW, Manna four-state and Zhang models, (a) 
i£[s|a] (average avalanche size for given avalanche area) vs. a, is presented, yielding ^ sa . (b) E[s\t] 
is given vs. t (avalanche time), yielding 7^. (c) i£[a|i] is given vs. t, yielding y at . System size is 
512 2 with 10 7 grains dropped. Data were binned with bin size increasing exponentially. 

FIG. 4. Typical large avalanche structure for the BTW model (a), Manna four-state model 
(b) and the Zhang model (c). Gray scales indicate the number of toppling events /(i) which 
occurred at each site during the avalanche. White represents zero relaxations, and black represents 
the maximal no. of relaxations [17 in (a), 34 in (b) and 22 in (c)]. System size is 128 2 . Note 
the shell structure in the BTW avalanche 0] vs. the irregular structure of the avalanche in the 
Manna four-state model, and the intermediate structure of the Zhang avalanche. These qualitative 
geometrical differences translate into quantitative differences in exponent values. 

FIG. 5. The activity profile of a typical large avalanche in the BTW, Manna four-state and 
Zhang model. The terrace number is plotted as a function of its area. The picture is reflected 
around the y axis. The system size is 128 2 . The BTW and Zhang models exhibit moderate slopes, 
while in the Manna four-state model the slope becomes extremely steep at high terrace numbers. 
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FIG. 6. The scaling functions for the BTW model for systemsize = 128, 256 and 512 (a) 
fu(v) which describes the time dependence of the energy during the avalanche, (b) fs(n) which 
describes the time dependence of the number of active sites and (c) which describes the time 

dependence of the avalanche area growth. We observe that all three scaling functions coincide, 
indicating a common scaling function for U, S and A. 

FIG. 7. The scaling functions for the Manna four-state model for systemsize = 128, 256 and 
512 (a) fuip)-, (b) /s(/-0 and (c) We observe that the scaling functions fu{f) and fs(v) 

coincide while /a(a*) has a completely different form. 

FIG. 8. The functions (a) fu{^), (b) /s(/x) and (c) for the Zhang model for 

systemsize = 128, 256 and 512. We observe that functions obtained from different system sizes do 
not coincide, indicating that these are not scaling functions. This indicates that the Zhang model 
lacks some of the characteristic features of a critical state found in the BTW and Manna models. 
Interestingly, for each system size fu(fJ-) and /s(/t) still coincide, while /a(/0 is different. 



19 




I I I I I I 1 1 I I I I I I 1 1 1 I I I I I I 1 1 1 I I I I I I 

BTW 
A Manna four-state 
Zhang 




+APo n 



A 



o. 



AAA A °o 
T +M O 

+ A co 



-i i i i i i ■ ■ ■ 



+ A co 
-HAAA coo 

jJ 1 1 — I 1 l ^fi'J'AO OOP 1 — 



io 2 io 3 io 4 io 5 

avalanche size s 



2.50 




2.00 



0.000 



0.005 



0.010 
1/L 



0.015 



0.020 



Fig. 2 



10' 



10 



10 



]*= 10 

io 2 



10 



10 







i i i 1 1 1 1 1 1 i i i i 1 1 1 1 1 i i i 1 1 1 1 1 1 i i i i 1 1 1 1 1 i i '^J^ 1 1 1 1 o 



I I I I I I I LI 



o BTW 

* Manna four-state 
+ Zhang 



+A 



A O 
s o 



iA O 

T O 

A o 



A o 
A O 
A ° 

^ o 



. A o 

A^ O 



At O 
A+ o° 

A+T© 



a£* 



-i i i i i i i i i i i i 



10° io 1 io 2 io 3 io 4 io 5 10* 



area a 



Fig. 3(a) 



10' 



10' 



& 10 



10 



10 







— I I I I I I I 1 1 I I I I I I I 1 1 I I I I I I 1 1 1 I I I I I I I 1 1 I I I I I I 

o BTW 

* Manna four-state ++ 
+ Zhang 



+ A* 

+ A 
. A 



4 



,ft 



ft 



ft 



ft 



ft' 



-i i i * i i i * i i i * i i i i i i i 



10' 



10' 



10 



10 



10 



10 



time t 



Fig. 3(b) 



10' 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I. 

• BTW 
JO 5 L a Manna four-state 
+ Zhang 



10 



— 3 

«, 10 

io 2 



10 



10 







** A 

4* * 



^ A 
* A 
A 



, * A 
lb A 
* A 



+ o 
* A 



. + ° A 
+ O A 

+ O A 

+ O A 

+ O A 



ft 



-i i i i i i i i i i i i i i i i i i ■ ■ 



10' 



10' 



10' 



10 



10 



time t 



Fig. 3(c) 



^^^^^^^^^^^ 




Fig. 4(a) 



(b) 




Fig. 4(b) 



(c) 




Fig. 4(c) 



35 



30 - 



25 - 



o BTW 

+ Manna four-state 
* Zhang 



+ 
+ 
+ 
+ 
+ 
+ 
+ 

-H- 
-H- 



O 

o 
H 



20 - 



15 - 



10 - 



5 - 





-15000 



hx xf 
■+o Of 
<* *> 

(S3- -6© 
«+ -H© 
<*+ +» 
«+ +® 
®+ +« 
30f -KX 

xa- Hex 

xe ex 
x-D <3-x 
x© ex 
xe ^x 
x© ex 

xe e< 

O 4x x+ O 
+ x x + 
i i i i i 



-5000 



5000 



15000 



Terrace size 



Fig. 5 




0.0 1.0 2.0 3.0 



t/(t) 



